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Abstract 

The 138 La (T 1/ / 2 =102 Gyr) - 138 Ce - 136 Ce system is proposed to be used as a nuclear cosmochronome- 
ter for measuring the time elapsed from a supernova neutrino process. This chronometer is applied 
to examine a sample affected by a single nucleosynthesis episode as presolar grains in primitive mete- 
orites. A feature of this chronometer is to evaluate the initial abundance ratio of 136 Ce/ 138 Ce using an 
empirical scaling law, which was found in the solar abundances. We calculate the age of the sample as 
a function of isotopic ratios, 136 Ce/ 138 Ce, and 138 La/ 138 Ce, and evaluate the age uncertainty due to 
theoretical and observational errors. It is concluded that this chronometer can work well for a sample 
with the abundance ratio of 138 La/ 138 Ce > 20 when the ratios of 136 Ce/ 138 Ce and 138 La/ 138 Ce are 
measured within the uncertainty of 20%. The availability of such samples becomes clear in recent 
studies of the presolar grains. We also discuss the effect of the nuclear structure to the v process 
origin of 138 La. 
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Long-lived radioactivities are used as nuclear cosmochronometers for dating the ages of nu- 
cleosynthesis episodes, metal-poor stars , and the Milky way |3j. Rutherford suggested 
the original idea, the uranium chronometer, in 1929 {4]. Recent progress of astronomical obser- 
vations has enabled to detect actinoid elements of early generations of stars [l, 2]. A fraction 
of the uranium in a metal-deficient star CS 31082-001 was measured and the stellar age was 
estimated with the 238 U chronometer for the first time in 2001 [1]. Only six chronometers with 
half-lives of 1-100 Gyr which are suitable for measuring cosmo logical ages were known. The 
other chronometers are 40 K and 232 Th suggested by Burbidge et al. [^j, 87 Rb and 187 Re sug- 
gested by Clayton p, and 176 Lu suggested by Audouze et al. It should be noted that these 
chronometers measure the age of slow or rapid neutron-capture reaction nucleosynthesis (s or r 
process). However, there is no cosmochronometer for neutrino- induced reactions in supernova 
explosions [y process). 
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The v process was suggested as the origin of heavy and light elements 



151 ]. and is of importance for studying neutrino spectra from the supernovae 



n 



ipj, 
11 
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m, 



and for 



discussing neutrino oscillation [15]. Among many heavy elements, only two isotopes of 138 La and 
180 Ta are considered to be synthesized primarily by the v process (8|, but there is another 
possibility for 180 Ta that is produced by the s process [if]]. The solar abundance of 138 La was 
quantitatively explained by the v process in theoretical calculations [ll, 17]. Thus, 138 La is a key 
for understanding the v process. In this paper, we propose a new cosmochronometer based on 
the abundance ratios of the three isotopes, 138 La, 138 Ce and 136 Ce, in a single sample affected 
strongly by a single nucleosynthesis episode. Recent studies of presolar grains in primitive 
meteorites provide samples affected strongly by single supernova nucleosynthesis episodes 18|, 
I19I ]. The isotopic patterns of heavy elements such as molybdenum and barium in the silicon 



carbide (SiC) presolar grains originating from the supernovae have been measured [20, Q. On 
the basis of these technological developments, the abundances of 136 > 138 Ce and 138 La in a single 
sample as the presolar grain are expected to be determined in future. 

The meta-stable isotope 138 La decays to 138 Ce or 138 Ba with a half-life of 102 Gyr. The 
138 La- 138 Ba (or 138 Ce) system has been used as a geochronometer to date the ages of rocks 221 ] . 
Audouze et al. pointed out that the 138 La system is not useful as the cosmochronometer since 
the abundance of the daughter uuc.eus popu.ated by the decay of -La may be miner relative to 
its initial abundance [7|. However, nucleosynthesis models cannot predict generally the initial 
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abundances of individual presolar grains, because the supernovae produce both of 138 La and 
138 Ce and their initial abundance ratios in individual grains are different. After the destroy 
of pre-existing 138 La and 138 Ce (and 136 Ce) by the weak s process in evolutionary stages of 



avers in the supernova explosion; 138 La 



progenitors, these nuclei are re-synthesized in different 
is synthesized by the v process in He- and C-rich layers 
7 process in O/Ne layers 23|, |24|. Meteorite compositions are produced in stellar outflows after 



ll| and 136 ' 138 Ce are synthesized by the 



mixture of materials originating from different layers, which may include both the Ce isotopes 
and 138 La. The mixing ratios of the different layers in the individual meteorite compositions are 
different. Therefore we should evaluate the initial abundances in each sample by an alternative 
method. 

Here we propose a novel empirical method to estimate the initial abundance of 138 Ce. A 
similar empirical method was suggested for a 146 Sm nuclear cosmochronometer, wherein an 
initial abundance of "*Sm is evaluated by the solar abundance pattern Q. In our proposed 
method, the initial abundance of 138 Ce can be calculated by an empirical scaling law for p- 
nuclei, which was found in the solar abundances 26J. The three isotopes, 138 La, 138 Ce, and 
136 Ce, are classified to the p-nuclei. The p-nuclei are located on the neutron-deficient side from 
the /9-stability line and their isotopic fractions are small (typically 0.1-1%). There are nine 
pairs of the two p-nuclei with the same atomic number. As shown in Fig. 1, the first and 
second p-nuclei are lighter than the s-nucleus by two and four neutrons, respectively. We found 
the empirical scalings that the abundance ratios of iV (lst p)/A^ (2nd p) and N Q (s)/N Q (lst 
p) are almost constant over a wide range of atomic number, respectively, where N & is the 
solar isotope abundance 26j. These scalings are a piece of evidence that the most probable 
origin of the p-nuclei is the supernova 7 process. In addition, we found a novel concept of "the 
universality of the 7 process", wherein the scalings hold for nuclides produced in individual 
supernova 7 processes under various astrophysical conditions [261] . We calculated the 7 process 
of core-collapse supernovae under various astrophysical conditions such as, metallicity, progen- 
itor mass, and explosion energy and found that the scalings hold for individual nucleosyntheses 
independent of the astrophysical conditions assumed 271 ] . We presented a detailed mechanism 
why the scalings fold for various supernovae in our previous papers 27, 2{J. Our proposed 
chronometer is based on the universality of the scaling for the two p-nuclei. 

Here we improve the empirical scaling for the two p-nuclei. There are nine pairs of the two 
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p-nuclei in the solar system, but we omit three elements of Mo, Ru and Er in Fig. [2] because 



the origin of p-nuclei of Mo and Ru is considered to be different from other p-nuclei [13 
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26] 



and a p-nucleus 164 Er is largely contributed from the s process 30j. The A^ Q (lst p)/A^ Q (2nd p) 
ratios in the other six elements increase with increasing the atomic number (see Fig. 2). We 
obtain a function, 

R PP (Z) = 0.0581 xZ - 2.18, (1) 

using a ^-square fitting. As shown in Fig. [2J the observed A^ (lst p) / N Q (2nd p) ratios are 
described well by this equation and the standard deviation of A^ Q (lst p) / N Q (2nd p) is only 9%. 

If abundance ratios of N( 138 La)/N( 138 Ce) and N( 136 Ce)/N( 138 Ce) in a sample affected 
strongly by a single supernova are known, the age T elapsed from the supernova episode is 
calculated by 

T = -^i 

ln2 

J Ney/N^c,) \ 

VAT( 1 38 La )/ A r(i38 Ce ) + Ih - J R pp (Ce)xAT( 136 Ce)/iV( 138 Ce)J / 

where N is the observed abundance in the sample, R pp is the initial AT (1st p)/N (2nd p) ratio 
at the freezeout of the 7 process, and b is the branching ratio of the j3~ decay of 138 La: \((3~) 
I (A(/3 _ )+A(EC)). The scaling holds for nuclides produced by individual supernova nucleosyn- 
theses and the function (JJJ would hold for individual presolar grains. The Rp P (Ce) ratio is 
estimated to be 1.19 using the function (JT]). Therefore, the age can be obtained from the two 
observed ratios of N( 138 La)/N( 138 Ce) and N( 136 Ce)/N( 138 Ce) in the sample. 

A part of 138 La may be contributed from the 7 process, but the age obtained from this 
chronometer is not affected by the contamination of the 7 process, because this chronome- 
ter measure the age elapsed from a supernova and both the 7- and v processes occur at the 
same supernova. The cosmic ray process [29j also has been suggested as the origin of 138 La. 
However, the presolar grains originating from the supernovae do not contain the comic ray 
products, because the pre-existing 138 La and 136 > 138 Ce originating from the interstellar media, 
which include the cosmic ray products, are destroyed by the weak s process before the super- 
novae. The p-nuclei that include the three isotopes are located on the neutron-deficient side 
from the /3-stability line, and hence they cannot be synthesized by neutron capture reactions in 
the weak s process. On the other hand, the pre-existing p-nuclei are destroyed by the neutron 
capture reactions. Therefore, the three isotopes in the presolar grains from the supernovae are 



newly synthesized by single nucleosyntheses without any pre-existing contamination such as 
the cosmic ray products. 

Here we discuss the uncertainty of the age using the chronometer. Figure 3 shows the cal- 
culated ages as a function of i? pp -N( 136 Ce)/N( 138 Ce) under the mixing ratio N( 138 La)/N( 138 Ce) 
= 1, 5, 10 and 20. We also calculated the ages as a function of N( 138 La)/N( 138 Ce) under the 
ratio N( 136 Ce)/N( 138 Ce) = 1.5, 2.0, and 3.0 (see Fig. 4). The uncertainty of the age is generally 
due to theoretical and observational errors. In the present proposed chronometer, the theo- 
retical errors come only from the uncertainty of the initial abundance ratio of R pp (Ce). The 
uncertainty of R pp (Ce) is estimated from the deviation of the observed iV (lst p)/iV Q (2nd p) 
ratios (see Fig. 2). The standard deviation of the observed ratios is 9% and we take 9% as the 
uncertainty of R pp (Ce). The observational errors are due to the errors of the N( 136 Ce)/N( 138 Ce) 
and N( 138 La)/N( 138 Ce) ratios in samples. Recently, measurement techniques with mass sepa- 
ration have been rapidly developed and the isotopic fractions of heavy elements such as Mo 
and Ba in SiC grains originating from the supernovae have been measured within the uncer- 



tainties of 7-15% 



20 



211 ] . Thus, we consider two cases that the observational errors of both 



the two ratios are 10 or 20%. Finally, we calculate the complete uncertainty due to both 
the theoretical and observational errors; namely, the errors of R pp , N( 136 Ce)/N( 138 Ce), and 
N( 138 La)/N( 138 Ce). If the ratios of N( 136 Ce)/N( 138 Ce) and N( 138 La)/N( 138 Ce) are measured 
with the uncertainty of 10%, respectively, the complete uncertainty is about 20%, 50% or 
90% in the case of N( 138 La)/N( 138 Ce) = 20, 10 or 5. With observational errors of 20%, the 
complete uncertainty turns out to be about 30% for a sample of N( 138 La)/N( 138 Ce) = 20. It 
should be noted that the complete uncertainty suddenly decreases with increasing the ratio of 
N( 138 La)/N( 138 Ce). 

Recent progress of the astronomical observation has enabled to detect actinoid elements 
of metal-poor stars. The ages of these metal-poor stars have been evaluated with U and Th 
chronometers within uncertainty of ~50% [l|, Q] • Therefore the present proposed chronometer 
is as useful as the cosmo chronometers have been used in the astronomy. 

The discussion mentioned above shows that the present proposed chronometer can work well 
if the presolar grains of N( 138 La)/N( 138 Ce) > 20 are found with the observational errors of 20%. 
A question we have to ask here is whether the presolar grains of N( 138 La)/N( 138 Ce) > 20 exist 
or not in nature. Since the solar abundance ratio A^ ( 138 La)/A^ ( 138 Ce) is about 0.14, it seems 
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to be difficult to find the samples of N( 138 La)/N( 138 Ce) > 20. However, recent studies of the 
presolar grains originating from the supernovae suggest that such samples are available. The 
carbon-rich grains of the SiC grains and low-density graphite grains are chemically condensed in 
a carbon-rich environment of C/O > 1. This suggests that the main component of the carbon- 
rich grains originating from the supernovae are produced in their C-rich layers, where 138 La is 
synthesized by the v process. In contrast, the contribution of O/Ne layers, where 136 ' 138 Ce are 
produced by the 7 process, is considered to be small. The isotopic fractions of several elements 
in individual presolar grains from supernovae were measured and theoretical calculations were 
performed to reproduce these isotopic fractions by proper mixture of materials produced in 



different layers |3l|, |32( . The observed isotopic fractions were reproduced consistently without 



the contribution from the O/Ne layers 3jj and with the small contribution from the O/Ne 



layers which is smaller than that from the C-rich layers by 2-4 orders of magnitude [32j. We 
stress again that 136 > 138 Ce are dominantly synthesized in the O/Ne layers by the 7 process. 
Therefore, we expect that samples with the abundance ratio of N( 138 La)/N( 138 Ce) > 20 will 
be found. Finally, we would like to stress the complete uncertainty suddenly decreases with 
increasing the ratio of N( 138 La)/N( 138 Ce). The discussion mentioned above shows that we can 
find samples, of which the N( 138 La)/N( 138 Ce) ratios are much larger than 20, and thereby the 
complete uncertainties of the ages of the sample become much lower than 30%. 

Although recent studies of the astrophysics suggest that 138 La may be dominantly synthe- 
sized by the v process the origin of 138 La has not been elucidated in astronomical 



observations or meteorite analyses. To take a hypothetical example, suppose that 138 La is 
dominantly synthesized by the cosmic ray process and the presolar grains are formed in an 
interstellar environment. The age estimated by this chronometer would be longer than that in 
the case that 138 La is dominantly synthesized by the v process because the cosmic ray nucle- 
osynthesis is a continuing process from the Galaxy formation to the solar system formation. 
Thus, the age measured in the presolar grains is of importance to understand the origin of 
138 La. 

Here we point out that the nuclear structure of 138 La may affect the performance of the 
chronometer and the v process origin of 138 La. A fact that the spin and parity of the ground 
state of an isotone 140 Pr (JV=81) is J 7r =l + suggests that a 1 + excited state may exist at a 



low excitation energy in 138 La. Although a 1 + state at 296 keV in 138 La was reported 



33|, 
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the lowest 1 + state has not been established experimentally. If this 1 + excited state exists at 
energy lower than 72 keV that is the energy of the first excited state, the 1 + state may be a (3 
unstable isomer since f3 decay can compete with an internal E4 transition. In such a case, our 
proposed clock cannot work well since the initial abundance cannot be estimated by using the 
scaling. In addition, this isomer affects the v process origin of 138 La because 138 La synthesized 
by the v process is destroyed via the isomer (see Fig. Ej). 

We would like to stress that the observed 1 + ground state of 140 Pr can not be understood 
by a simple shell- model picture that would lead to the ground state with J n =2 + or 3 + 34J. 
This indicates that there is a mechanism decreasing the energy of the 1 + state. In order to 
examine the existence of the 1 + isomer, we calculate the nuclear structures of 138 La and 140 Pr 
in a shell-model. In the mass region around 138 La, the last proton (neutron) can occupy either 
d.5/2 or g7/2 ^3/2) orbit around the Fermi surface. The 1 + and 5 + states in 138 La and 140 Pr can 
be understood generally by (vrd 5 / 2 )®(^d3/2) _1 and (vrg 7 / 2 )®(i/d3/2) _1 configurations, respec- 
tively, where n (v) is proton (neutron). Figure [6] show the calculated results. We successfully 

351 ] and 140 Pr 36J lower than 500 keV and, in 



34, 



describe the observed excited states in 138 La 
particular, the ground state spin of 140 Pr. We find that the coupling of two- and four-particle 
excited configurations of 7r(g7/ 2 )~ 2 (d 5 /2) 3 and vr(g 7 /2)~ 4 (d 5 /2) 5 leads to the 1 + ground state. A 
calculated 1 + state in 138 La locates at 169 keV. This result suggests that the 1 + state in 138 La 
is not an isomer but an experimental determination is desired. 

In summary, we propose a novel 138 La (Ti/2=102 Gyr) - 138 Ce - 136 Ce nuclear cos- 
mochronometer for measuring the time elapsed from the supernova v process. We introduce 
a novel method using the universality of the empirical scaling for p nuclei to evaluate the 
initial abundance of 138 Ce. The chronometer is applied to examine a sample affected by single 
supernova nucleosyntheses as presolar grains in primitive meteorites. Recent studies of the 
presolar grains provide isotopic fractions of several different elements of single samples. The 
chronometer can work well for a sample, of which its abundance ratio of 138 La/ 138 Ce is larger 
than 20, when the ratios of 136 Ce/ 138 Ce and 138 La/ 138 Ce are measured within the uncertainty 
of 20%. The availability of such samples becomes clear in recent studies of the presolar 
grains. With the chronometer, we can measure both the age of the supernova associated 
with neutrino-driven wind and the isotopic abundance patterns of elements synthesized by 
the supernova. We present that the energy of the lowest 1 + state may affect the chronometer 
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performance and the v process origin of 138 La, and its measurement is desired. 
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FIG. 1: A partial nuclear chart and nucleosynthesis flows around 138 La. The nucleus 138 La is shielded 
by stable isobars against (3 decay and is produced by neither s nor r process. In the v process, 138 La 
is mainly synthesized from 138 Ba via the charged current reaction. Two p-nuclei 136 Ce and 138 Ce are 
mainly produced from an s-nucleus 140 Ce via successive (7,n) reactions. 




55 60 
Atomic Number 



70 



FIG. 2: Observed isotope ratios of the first p- nuclei to the second p- nuclei with the same atomic 
number. The solid line is the result of x-square fitting. The standard deviation of the fitting is 9%. 
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FIG. 3: Calculated age as a function -R pp -N( 136 Ce)/N( 138 Ce). The age is the time elapsed from a 
nucleosynthesis episode to the present, which includes the solar system age 4.6 Gyr. The solid curves 
of (a), (b), (c) and (d) are the calculated ages with mixing ratios N( 138 La)/N( 138 Ce) = 20, 10, 5 and 
1, respectively. The dot-dashed line is the age of T=10 Gyr. The solid, dashed, and dot lines are 
the upper and lower limits arisen from uncertainty of R pp or N( 136 Ce)/N( 138 Ce) for (a), (b) and (c), 
respectively. 




FIG. 4: Calculated age as a function N( 138 La)/N( 138 Ce). The age is the time elapsed from a 
nucleosynthesis episode. The solid curves of (a), (b) and (c) are the calculated ages with the ra- 
tios N( 138 Ce)/N( 136 Ce) = 1.5, 2.0 and 3.0, respectively. The dot-dashed line is the age of T=10 
Gyr. The solid, dashed, and dot lines are the upper and lower limits of the assumed uncertainty of 
N( 138 La)/N( 138 Ce) for (a), (b) and (c), respectively. 
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Low s P in Excited 
I states 




138 Ba 

FIG. 5: Synthesis and destruction of 138 La in the v process. Low spin excited states such as J n =l + 
are strongly populated from 138 Ba by Gamow- Teller transitions. If the 1 + isomer exists, these excited 
states strongly decay to the 1 + isomer which decays to daughter nuclei. 
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FIG. 6: The calculated and measured levels in 140 Pr and 138 La. The observed levels of 140 Pr and 
138 La are taken from ref. and refs. 34, 



351 ]. respectively. 
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